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Abelian mechanism of non-Abelian color confinement is observed in a gauge-independent way 
by high precision lattice Monte Carlo simulations in gluodynamics. An Abelian gauge field is 
extracted with no gauge-fixing. A static quark-antiquark potential derived from Abelian Polyakov 
loop correlators gives us the same string tension as the non-Abelian one. The Hodge decomposition 
of the Abelian Polyakov loop correlator to the regular photon and the singular monopole parts also 
reveals that only the monopole part is responsible for the string tension. The investigation of the flux- 
tube profile then shows that Abelian electric fields defined in an arbitrary color direction are squeezed 
by monopole supercurrents with the same color direction, and the quantitative features of flux 
squeezing are consistent with those observed previously after Abelian projections with gauge fixing. 
Gauge independence of Abelian and monopole dominance strongly supports that the mechanism of 
non-Abelian color confinement is due to the Abelian dual Meissner effect. 
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Color confinement in quantum chromodynamics 
(QCD) is still an important unsolved problem 
't Hooft and Mandelstam 0] conjectured that the 
QCD vacuum is a kind of a magnetic superconducting 
state caused by condensation of magnetic monopoles and 
an effect dual to the Meissner effect works to confine 
color charges. However, in contrast to SUSY QCD Q 
or Georgi-Glashow model 0, Q with scalar fields, to find 
color magnetic monopoles which condense is not straight- 
forward in QCD. 

An interesting idea to realize this conjecture is to 
project SU(3) QCD to an Abelian [U(l)] 2 theory by a 
partial gauge fixing Q. Then color magnetic monopoles 
appear as a topological object. Condensation of the 
monopoles causes the dual Meissner effect [tl fic| . 
However there are infinite ways of the above partial 
gauge-fixing and whether the 't Hooft scheme is gauge 
independent or not is not clear. Moreover why non- 
Abelian color charges are confined in the framework of 
the Abelian mechanism is not clarified. 

Numerically, an Abelian projection in non-local gauges 
such as the maximally Abelian (MA) gauge [ll|, Q3, [l3| 
has been found to support the Abelian confinement sce- 
nario beautifully [3, M, M, EE M, 03 HI- Very re- 
cently, the present authors have shown that the Abelian 
dominance and the dual Meissner effect are observed 
clearly also in local unitary gauges such as F12 and 
Polyakov (PL) gauges [2l| . These results strongly sug- 
gest that the Abelian confinement mechanism is gauge- 
independent. 

In this Letter, we study the QCD vacuum after ex- 
tracting an Abelian link field in a completely gauge- 
independent way without adopting any special local or 
non-local gauge fixing. We observe that an Abelian con- 
finement mechanism due to condensation of monopoles 
is realized. A static potential derived from Abelian 
Polyakov loop correlators gives us the correct string ten- 



sion. Moreover only the monopole part in the Abelian 
Polyakov loop is responsible for the string tension. 
Abelian electric fields defined in an arbitrary color di- 
rection are squeezed and the corresponding monopole 
currents play the role of magnetic supercurrents. States 
which are neutral in all color directions are not confined 
and appear as a physical state. It is just a color-singlet 
state. Hence, confinement of non-Abelian color charges, 
not that of Abelian charges, is explained in the frame- 
work of the gauge-independent Abelian mechanism. 

These findings are completely novel and exciting, al- 
though the continuum and the infinite-volume limits are 
not studied. 
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FIG. 1: The Abelian static potential in comparison with the 
non-Abelian one. The solid lines denote the best fit to a 
function Vm- 

Firstly we discuss an Abelian static potential. We gen- 
erate thermalized gluon configurations using the Wilson 
action at a coupling constant (3 = 2.5 on the lattice A^ 4 — 
24 4 , where the lattice spacing a(/3) = .0836(8) [fm]. For 
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simplicity we consider SU(2) gluodynamics, since essen- 
tial features are not altered in SU(3). We extract a 2 x 2 
diagonal Abelian link field in an arbitrary color direction. 
For example, in the 173 direction, 



M s ) = U°Js) + iaU^s) = C M (s)-diag 
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where # M (s) = arctan(J7^(s)/J7°(s)). Note that we can 
do the same also in the a\ or 172 direction, since all three 
components are equivalent with no gauge-fixing. By us- 
ing the multi- level noise reduction method [22| , we evalu- 
ate the Abelian static potential from the correlation func- 
tion of the Abelian Polyakov loop operator 
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separated at a distance R. For the multi-level method, 
the number of sublattices adopted is 6 and the sublat- 
tice size is 4. The results are surprisingly beautiful as 
seen from Fig. [TJ To reduce the lattice artifact due to 
finite-lattice cutoff, we plot the potential using 0{a 2 ) 
improved distances [HI, [24| ■ We try to fit the data to 
a usual function Vgt = crR — c/R + fj, and find almost 
the same string tension and the Coulombic coefficient as 
shown in Table HI indicating Abelian dominance. Here 
the number of independent vacuum configurations is 10 
in all cases. The errors are determined by the jackknife 
method. Our results of the string tension are consistent 
with theoretical observations on the basis of reasonable 
assumptions [25l.[2^|. 
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FIG. 2: The static potential (top) and the force (bottom) 
from the Abelian (Pa), the monopole (Pmon) and the photon 
contributions (P p h) in comparison with the non- Abelian ones. 



by a product of a photon and a Dirac-string parts [27 
Note that 



TABLE I: Best fitted values of the string tension a 2 a, the 
Coulombic coefficient c and the constant a/i. NA and A-NGF 
denote Non-Abelian and Abelian with no gauge-fixing. Mup 
is the number of internal updates in the multi- level method. 
FR means the fitting range. The x 2 for the central value is 
X 2 /N df < 0.1. 





aa 2 


c 




FR(R/a) 


Mup 


NA 
A-NGF 


0.0348(7) 
0.0352(16) 


0.243(6) 
0.231(39) 


0.607(4) 
1.357(17) 


3.92 - 9.97 
4.94 - 9.97 


15000 
160000 



9 4 (s) = -Y,D( S -s')[dX 4 (s') + d 4 (dl9„(s'))} , (2) 

s' 

where D(s — s') is the lattice Coulomb propagator, 
0,j,u(s) = dfj,6 v (s) - dyO^s) and d v (d' v ) is a for- 
ward(backward) difference. We have used d v d' v D(s — 
s 1 ) = —§ ss i. The second term in the right-hand side 
of ([2]) does not contribute to the Abelian Polyakov loop 
Now 6> M „(s) = 9^(s) + 27rrv,(s) (|0 M „| < tt), where 
n lJlv (s) is an integer corresponding to the number of the 
Dirac string. Hence we obtain Pa = P p h ■ Pmon, where 



TABLE II: Best fitted values of the string tension aa 2 and 
the Coulomb coefficient c. M-NGF (P-NGF) denotes the 
monopole (the photon) part. 





aa 2 


c 


/ia 


FR(R/a) 


X 2 /N df 


NA 


0.181(8) 


0.25(15) 


0.54(7) 


3.92 - 8.50 


1.00 


A-NGF 


0.183(8) 


0.20(15) 


0.98(7) 


3.92 - 8.23 


1.00 


M-NGF 


0.183(6) 


0.25(11) 


1.31(5) 


3.92 - 6.71 


0.98 


P-NGF 


-0.0002(1) 


0.010(1) 


0.48(1) 


4.94 - 9.44 


1.02 



Secondly we discuss the role of monopole contribution. 
The monopole part of the operator can be extracted as 
follows. The Abelian Polyakov loop §T§ can be written 



N-l 

P ph = exp{^^^T;( S + fc4- S XMs')}, 

fc=0 s' 
N-l 

Pmon = exp{-27ri^^D(s + fc4- s')d' v n v4 ,{s')}. 

k=0 s' 

We call Pph and P mon the photon and the monopole 
contributions, respectively, since the Dirac string 
n,j3 1 (s) leads us to a monopole current fc M (s) = 
(l/2)e /1Q , J 8 7 9 a n / 9 7 (s + /i) 

We need a non-local Coulomb propagator in the sep- 
aration, so that the multi-level noise reduction method 
cannot be applied in this case. Here we consider aT^O 
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system in the confinement phase with the Wilson ac- 
tion on 24 3 x 4 lattice. We use about 6000 thermal- 
ized configurations at (3 = 2.2, where the lattice spacing 
is a(f}) = .191(8) [fm]. Since the expectation values of 
the correlation functions of Pa, P p h and Pmon are still 
very small with no gauge-fixing, we adopt a new noise 
reduction method. For a thermalized vacuum ensemble, 
we produce many gauge copies applying random gauge 
transformations, compute the operator for each copy, and 
take the average of all copies. Note that as long as a 
gauge-invariant operator is evaluated, such copies are 
identical, but they are not if a gauge-variant operator 
is evaluated. Practically, we prepare 1000 gauge copies 
for each configuration. We also apply one-step hypercu- 
bic blocking (HYP) [29| to the temporal links for further 
noise reduction. 

We obtain very good signals for the Abelian, the 
monopole and the photon contributions to the static po- 
tential as shown in Fig. O We try to fit the poten- 
tial in Fig. [2] to the function Vat and extract the string 
tension and the Coulombic coefficient of each potential 
as summarized in Table [TTJ Abelian dominance is seen 
again beautifully in this case. Moreover, we can see 
monopole dominance, namely, only the monopole part of 
the Polyakov loop correlator is responsible for the string 
tension. The photon part has no linear potential. The 
agreement among the string tensions coming from non- 
Abelian, Abelian and monopole Polyakov loop correla- 
tors is almost perfect as seen also from the force in Fig. [5] 
in comparison with the MA case, where only 80-90 per- 
cent agreement is observed at finite lattice spacings. The 
short-range behavior of the potential may be affected by 
HYP. 

Thirdly we discuss the Abelian dual Meissner effect. 
We investigate the Abelian flux-tube profile by evaluat- 
ing connected correlation functions [3(| HH between a 
Wilson loop W and Abelian operators Oa constructed 
by Abelian link fields, 



(0 A (r)}w = 



(Tr [LW{r = 0, R, T)L^a 3 A {r)\ ) 
(Tr \W(R, T)]) 



where L is a product of non-Abelian link fields (a 
Schwinger line) connecting the Wilson loop with the 
Abelian operator. We may use the cylindrical coordi- 
nate (r, 4>, z) to parametrize the the q-q system, where 
the z axis corresponds to the q-q axis and r to the trans- 
verse distance. We are interested in the field profile as a 
function of r on the mid- plane of the q-q distance. In this 
calculation, we employ the improved Iwasaki gauge ac- 
tion [32| with the coupling constant f3 = 1.20, which cor- 
responds to the lattice spacing a(J3) = .0792(2) [fm] [33| . 
The lattice volume is 32 with periodic boundary condi- 
tions. We generate 4000 thermalized configurations. To 
improve a signal-to-noise ratio, the APE smearing tech- 
nique is applied to the Wilson loop [34j . 

We measure all components of the Abelian electric 
fields Eaz{s) = 04i(s) and find that only Ea z is squeezed 
as shown in Fig. [31 We try to fit (Ea z )w to a function 
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FIG. 3: The profile of the Abelian electric fields for W(R 
5a,T = 5a). 
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FIG. 4: The curl of the Abelian electric field, magnetic 
displacement currents and monopole currents for W(R = 
5a,T = 5a). 

f(r) = ci exp(— r/X) + c$. Here A corresponds to the 
penetration length. We obtain A = 0.128(2) [fm], which 
is similar to those obtained in the MA gauge and unitary 
gauges [2l| as seen from Table ??. 
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FIG. 5: The correlation between the Wilson loop and the 
squared monopole density for W(R = 5a, T — 5a). The solid 
line denotes the best exponential fit. 

To see what squeezes the Abelian electric field, let us 
study the Abelian (dual) Ampere law 

V x E A = d 4 B A + 27rfc , 

where BAi(s) — (l/2)e i ji,()ji i .(s). Each term is evaluated 
on the same mid-plane as for the electric field. We find 
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that only the azimuthal components are non- vanishing, 
which arc plotted in Fig. 2J Note that if the electric field 
is purely of the Coulomb type, the curl of electric field 
is zero. Contrary, the curl of the electric field is non- 
vanishing and is reproduced only by monopole currents. 
The magnetic displacement current is almost vanishing. 
These behaviors are clearly a signal of the Abelian dual 
Meissner effect, which are quit e the same as those ob- 
served in the MA gauge [Ulloj]. 

Fourthly, we may estimate the vacuum type by eval- 
uating also the coherence length £ from the correla- 
tion function between the Wilson loop and the squared 
monopole density k^(s) [35| . The correlation function 
is plotted in Fig. \E\ and the coherence length extracted 
from the functional form g(r) = c[ exp(— y/2r/£) + c' 
is £/V2 = 0.102(3) [fm]. The GL parameter V2k = 
A/£ = 1.25(6) is close to the values obtained with gauge 
fixing (2l|. Since the Wilson loop used here may still be 
small, what we can say is that the vacuum type is near 
the border between the type 1 and 2. 

To summarize, we have observed gauge-independence 
of Abelian and monopole dominance for the string ten- 
sion and of the Abelian dual Meissner effect in gluody- 
namics by using lattice Monte Carlo simulations. These 
results are quite remarkable in the sense that confine- 
ment of non- Abelian color charges can be explained in 
the framework of the Abelian dual Meissner effect. Since 
no gauge- fixing is done, gauge fields in any color direc- 



tion are equivalent. Abelian electric fields in all color 
directions are squeezed due to monopoles in the corre- 
sponding color direction. An Abelian neutral state in 
all color directions can survive as a physical state, and 
such a state is only the color singlet state. For exam- 
ple, consider meson states u c u c and d c d c , where u c (d c ) 
is an eigenstate of CT3/2 with an eigenvalue 1/2 (—1/2). 
These are Abelian neutral in the 03/2 direction. Sim- 
ilarly, U C U C and D C D C are Abelian neutral in the cri/2 
direction, where U c = (u c + d c )/ \]~2 (D c — (u c — d c ) / %/2) 
is an eigenstate of eri/2. Note that u c u c (U C U C ) and d c d c 
(D C D C ) contain both Abelian charged and neutral states 
in the <J\/2 (03/2) direction. But a SU(2) singlet state 
u c u c + d c d c = U C U C + D C D C is Abelian neutral in all 
color directions. Hence confinement of non- Abelian color 
charge can be explained in terms of the Abelian confine- 
ment scenario of the dual Meissner effect. 

Finally it is interesting to study the relation between 
the violation of the non- Abelian Bianchi identity (36| and 
Abelian monopoles with no gauge-fixing. 
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